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Volume dependence of light hadron masses in full lattice QCD 

B. Orth a , T. Lippert a and K. Schilling 3 - 
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The aim of the GRAL project is to simulate full QCD with standard Wilson fermions at light quark masses 
on small to medium-sized lattices and to obtain infinite-volume results by extrapolation. In order to establish 
the functional form of the volume dependence we study systematically the finite-size effects in the light hadron 
spectrum. We give an update on the status of the GRAL project and show that our simulation data for the light 
hadron masses depend exponentially on the lattice size. 
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1. INTRODUCTION 

Lattice QCD has always been plagued by the 
fact that towards realistic quark masses the sim- 
ulation costs increase drastically due to the large 
correlation lengths of the light states and the 
large lattice volumes that one employs in order 
to avoid finite-size (FS) effects. This difficulty 
has, in fact, prevented us until this day from car- 
rying out simulations with truly physical param- 
eters. Instead one resorts to computations with 
several un-physically large quark masses in FS- 
effect-free volumes, followed by an extrapolation 
to the physical quark mass. However, the chiral 
extrapolation introduces uncertainties that can 
be reduced only by increasing the overlap of lat- 
tice simulations with ChPT. It is thus inevitable 
to look out for feasible ways to obtain lattice re- 
sults at light quark masses. Current approaches 
like overlap, domain-wall or improved staggered 
fermions are still very costly, especially if one goes 
beyond the quenched approximation. Faced with 
limited computer resources we therefore investi- 
gate the feasibility of an approach that restricts 
itself to the use of standard Wilson fermions. 

The aim of the GRAL project is to carry out 
a systematic study of FS effects in the light 
hadron spectrum in the parameter regime acces- 
sible to us. With simulations on series of small to 
medium-sized lattice volumes (at fixed coupling 
and quark mass, respectively) we also want to 
address the issue of whether one can use formu- 
lae obtained from this investigation to extrapo- 
late hadronic observables to the infinite volume. 



2. VOLUME DEPENDENCE OF THE 
LIGHT HADRON MASSES 

In 1983 M. Liischer published a universal field 
theoretic formula for the shift in the mass of a sta- 
ble particle enclosed in a finite box |1I2| . The for- 
mula states for asymptotically large volumes that 
the FS mass-shift vanishes exponentially with in- 
creasing box size at a rate that depends on the 
particle considered and on the spectrum of light 
particles in the theory. If, in QCD, we consider 
a stable hadron in a box of size L 3 with periodic 
boundary conditions and sufficiently large (Eu- 
clidian) time-extent T, its mass Mh{L) becomes 
a universal function of M n L (where, at fixed (3, 
M n = linif^oo M n (L)) in the finite- volume con- 
tinuum limit. Since FS effects probe the system 
at large distances and are thus insensitive to the 
form and magnitude of the UV cutoff, this func- 
tion is expected to hold also for finite lattice spac- 
ings. With the (constant) LO chiral expression 
for the 7r-7r forward scattering amplitude the rel- 
ative mass-shift for the pion is given by 
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If we assume, for simplicity, the pion-nucleon 
forward scattering amplitude to be constant, 
too, then the formula for the nucleon involves 
a term that is — up to low-energy constants — 
proportional to the one for the pion. In addition 
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it contains a term proportional to 

exp (-M„Ly/l - M^/(AM%)j /(M„L). (3) 

On the other hand, in 1992 Fukugita et al. 
found a power-like L-dependence of the form 

M H (L)=M H + c/L 3 (4) 

in their data for the pion, rho and nucleon 
masses |3I4| . The usual explanation for the ap- 
parent discrepancy is that Liischer's formula deals 
with asymptotically large lattice volumes where 
FS effects arise from a squeezing of the vir- 
tual pion cloud surrounding the hadron in a box 
with periodic boundary conditions. In contrast, 
Fukugita et al. ascribe the power-like behaviour 
of their data to a squeezing of the hadron itself 
(as it would occur for small lattice volumes). In 
this picture one expects a power-like dependence 
for sub-asymptotic lattice volumes that gradually 
changes into an exponential fall-off of the hadron 
masses towards larger box sizes. 

3. GRAL STATUS 

Table ^ summarises some details of our simu- 
lations. We use the standard unimproved Wil- 
son action with two degenerate flavours of dy- 
namical fcrmions. Our HMC codes are running 
on the APEmille at DESY/Zcuthcn and on the 
cluster installation ALiCE at Wuppertal Univer- 
sity. In addition to the already available SESAM 
(L = 16) and T%L (L = 24) ensembles we 
are generating configurations on smaller volumes 
with L varying between 10 and 16. The figures 
referring to these runs are, where available, pre- 
liminary. With our simulations at (3 — 5.32, 
k = 0.1665 we aim at a quark mass correspond- 
ing to Mir/Mp « 0.5, i.e. below the smallest 
M n /M p of 0.575(16) previously achieved by the 
SESAM/T^L project. In this article we focus on 
the ensembles at (3 = 5.6, n = 0.1575 where we 
assume the largest, 24 3 x 40-lattice to be free of 
FS effects. 

4. RESULTS 

The (preliminary) masses we obtain for the 
pion, rho and nucleon are displayed in Fig. ^ 



1 


1 

t 


I | i | i 


0.8 






0.6 


t 


• N 


0.4 




P 




0.2 


1,1,1,1, 



I i I i I i I i L_ 

8 12 16 20 24 



Figure 1. Dependence of the n-, p- and A-masses 
for (3 = 5.6, k — 0.1575 on the spatial lattice 
extent, L. The curves are explained in the text. 



The flat dashed curve represents Liischer's for- 
mula for the pion mass, Eq. QJ, with M w and 
as obtained from the 24 3 -lattice. From this 
curve it is obvious that the data points at L = 
10, 12, 14 are outside the chiral regime. This can 
also be seen from F^L, which takes the values 

0. 91(3), 1.10(3), 1.28(4), 1.46(4) and 2.19(6) for 
L = 10, 12, 14, 16, 24, respectively. Since ChPT 
is a low-energy (large distance) expansion it is 
only valid for F^L > 1 and M^/^AnF^) <C 1; 
M„/(4ttF„) is 0.240(7) for L = 24. Nevertheless 
we observe that the pion data can be well fitted 
from L — 12 with the "semi-empirical" formula 

M m m 3 ml' 2 e~ a ^ L 
M 7r (L) = m n + a 1 ^-^ —^ (5) 

where m„, a\ and a-i are free fit-parameters. a\ 
and d2 are chosen such that their deviation from 1 
indicates "how far away" we are from the original 
formula 0J. For the rho (which cannot decay on 
the employed lattices) and the nucleon we adopt, 
for various practical reasons, the generic ansatz 

g — a 2 m 7r L 

M PtN (L) = m PtN + ai L3/2 , (6) 

1. e. we effectively drop the term proportional 
to Here, is fixed to the value obtained 
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Table 1 

GRAL production status as of 07/08/2003 






K 


L 3 x T 


A con f 


<±TrD> 




o" 1 [GeV] 


M n (L)L 




L [fm] 


5.32 


0.1665 


12 3 x 32 
14 3 x 32 


4900 
2300 


0.5395(4) 
0.5382(2) 
















16 3 x 32 


3400 


0.5382(1) 


3.84(3) 


1.51(2) 


4.3(2) 






5.6 


0.1575 


10 3 x 32 


6200 


0.57324(7) 






4.9(2) 


2.76(3) 


0.850(4) 






12 3 x 32 


7200 


0.57280(5) 






4.41(8) 


3.32(4) 


1.020(5) 






14 3 x 32 


6500 


0.57262(4) 






4.23(7) 


3.87(4) 


1.190(6) 




SESAM 


16 3 x 32 


6500 


0.57254(3) 


5.96(8) 


2.35(3) 


4.49(6) 


4.43(5) 


1.360(7) 




TxL 


24 3 x 40 


5100 


0.57248(1) 


5.89(3) 


2.32(1) 


6.64(7) 


6.64(7) 


2.04(1) 




0.1580 


14 3 x 32 
16 3 x 32 


4400 
3900 


0.57367(4) 
0.57345(5) 






4.0(2) 
3.7(2) 








TxL 


24 3 x 40 


4500 


0.57337(2) 


6.23(6) 


2.45(2) 


4.78(7) 







from the pion fit. All the fits are displayed in 
Fig. 2] where the lines are solid within the fit 
range and dotted outside. For comparison the 
dash-dotted curves show 1 / L 3 -fits of the form (JIJ . 
The x 2 /dof- values for the exponential fits are 
0.20, 0.03 and 1.33 for pion, rho and nucleon, re- 
spectively. In contrast, for the 1/V-fits x 2 /dof is 
19, 8 and 4. The asymptotic mass (m p , rajy) 
obtained from the fit deviates less than 1% from 
M„ (M p , M N ) at L = 24 and lies thus well within 
the errors of M n (M p , M N ). 

5. CONCLUSIONS AND OUTLOOK 

We clearly see an exponential behaviour in our 
data for the light hadrons, as opposed to the 
l/L 3 -dependence reported earlier by Fukugita 
et at. However, the l/L 3 -fits are acceptable 
when we shift the fit-range from (12. .24) towards 
smaller L, say (10. .16). Moreover, we have to 
allow for extra free parameters in Liischer's for- 
mula in order to achieve a good description of our 
data. These parameters are 0(1000) in the case 
of ai and about 2 in the case of 02 (except for 
the nucleon where a\ = O(100) and 0,2 ~ 1) and 
vary with the left fit boundary. This, together 
with the F^L-values mentioned above, suggests 
that we are in an intermediate, sub-asymptotic 
regime were we merely see the onset of the chiral 
behaviour predicted by Liischer's formula. 



In a next step we will analyse our data at 
= 5.6, k = 0.158 and = 5.32, k = 0.1665 
in order to check whether the observed volume- 
dependence is reproduced there. We will then 
address the question if and to what precision an 
extrapolation to the infinite volume is possible. 
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